Signatures of fractal clustering of aerosols advected under gravity 
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Aerosols under chaotic advection often approach a strange attractor. They move chaotically on 
this fractal set but, in the presence of gravity, they have a net vertical motion downwards. In practical 
situations, observational data may be available only at a given level, for example at the ground level. 
We uncover two fractal signatures of chaotic advection of aerosols under the action of gravity. Each 
one enables the computation of the fractal dimension Do of the strange attractor governing the 
advection dynamics from data obtained solely at a given level. We illustrate our theoretical findings 
with a numerical experiment and discuss their possible relevance to meteorology. 
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The transport of finite-size particles plays important 
role in several fields, from cloud physics [1] to plankton 
dynamics [2]. The recent interest [3, 4] in this problem 
comes in part from the fact that the dynamics of these 
particles is dissipative due to the drag force. This makes 
dynamical systems tools and concepts, such as attractors 
and dimensions, applicable. In the presence of gravity, a 
net vertical motion occurs due to the density difference 
between fluid and particle. For heavy particles [aerosols), 
this leads to raindrop falling in the atmosphere [5] and 
to the sedimentation of plankton [2] and marine snow [6] 
in the ocean. 

In such situations, knowledge of the advection dynam- 
ics of the aerosols is of fundamental importance, whereas 
usually only data obtained at a given level (height) are 
available. This occurs often, for instance, in meteorology, 
in the case where the aerosols are raindrops. In fact, it 
is much easier to obtain direct measurements of the rain- 
drops when they reach the ground level than before, i.e., 
when they are being advected in the air flow. The deriva- 
tion of approaches to obtain information on the advection 
dynamics of aerosols solely from data observed at a given 
level is, therefore, an instrumental and relevant task. In 
particular, here we are interested in approaches to obtain 
the fractal dimension of the set where the aerosols cluster 
while they are advected. 

We report the uncovering of two independent frac- 
tal signatures of chaotic advection under gravity. Both 
make the computation of the fractal dimension of 
the strange attractor in the A^-dimensional configuration 
space possible without prior knowledge of the advection 
dynamics. First, we show that the time series of the in- 
stants of arrival of advected aerosols in a small detector 
placed at a given level has a fractal dimension which is 
equal to 

do^\ + Do-N. (1) 

We assume that Dq < N, which implies that the at- 
tractor in the full 2iV-dimensioiial phase space is Dq- 
dimensional since a set of dimension Dq, when pro- 
jected into a space of dimension N, typically remains 



Do-dimensional ii < N [7]. Second, we show that 
the spatial distribution of the aerosols reaching a line at 
a given level contains discontinuities (jumps) at points 
that form a fractal set whose dimension is again equal 
to c?o [8]- We illustrate our findings with a numerical 
experiment. 

The dimensionless form of the governing equation for 
the path r(t) of aerosols much denser than the fluid, sub- 
jected to Stokes drag and gravity, reads as [9]: 

r = A (u - r - Wn) , (2) 

where r is the velocity of the aerosol, u = u(r(t), t) is the 
fluid velocity fleld evaluated at the position r{t) of the 
aerosol, and n is a unit vector pointing upwards in the 
vertical direction. Throughout this paper we consider the 
vertical direction along the axis y. The inertia parame- 
ter A (larger values for smaller inertia) can be written 
in terms of the densities pp and p/ of the aerosol and 
of the fluid, respectively, the radius a of the aerosols, the 
fluid's kinematic viscosity v, and the characteristic length 
L and velocity U of the flow. It is A = R/ St, where 
R = Pf/pp <^ 1 and St = {2a^U)/{%vL) is the Stokes's 
number of the aerosol. As seen from (2), the gravita- 
tional parameter W provides the dimensionless settling 
velocity in a medium at rest. The actual settling velocity 
is the result of two effects: the gravitational attraction 
(buoyancy) and an updraft, if present. We consider the 
settling velocity to be comparable with U , implying a W 
of the order of unity. 

For convenience, we treat the case where the fluid 
flow is two-dimensional, iV = 2. In this situation, the 
phase space of the advection dynamics of the aerosols is 
4-dimensional, since the aerosols are not constrained to 
move with the same velocity as their corresponding fluid 
elements. For the sake of concreteness, let us consider 
the time-smoothened version of the alternating sinusoidal 
shear flow of Ref. [10]. In dimensionless form it is given 
as 

f Ux{Y,t) = 0.5(1 -|-tanh(7sin(27ri)))sin(27ry), 

(3) 

{ Uy{r,t) = 0.5(1 - tanh(7sin(27rt)))sin(27r.T). 
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FIG. 1: Snapshots, taken at (a) t = 0.65, (b) 0.75, (c) 0.85 
and (d) 0.95, of the projection A of the strange attractor into 
the configuration space. 



This flow is defined on the unit square with periodic 
boundary conditions, and is periodic in time with a unit 
period. The vertical direction corresponds to the y-axis. 
Apart from the spatial (sinusoidal) factor, each veloc- 
ity component consists of two plateaus in time with a 
rapid but smooth crossover if 7 = 20 /tt. This is a simple 
analytically given model which, nevertheless, possesses 
a paradigmatic property of both chaotic and turbulent 
flows [11]: the intense stretching of material elements. 

Substituting Eq. (3) into Eq. (2) and fixing R = 
10^^ and W = 0.8, there are regions of the parameter 
St for which the dynamics of the aerosols is ruled by a 
strange attractor. We note that the timc-indcpcndence 
of the attractor's dimension Dq, which will be essential in 
what follows, is valid for a broad class of randomly time- 
dependent flows as well (see conclusion). A bifurcation 
diagram and our analysis suggest the existence of strange 
attractors of dimension between 1 and 2 in an interval 
around St = 2 x 10~^ (apart, as usual, from some small 
periodic windows). Here wc illustrate our general ideas 
by analyzing the strange attractor corresponding to St ~ 
2 X 10"'*, which is prototypical. Figure 1 shows four 
snapshots of the projection. A, of the strange attractor 
into the configuration space. Because of gravity, there 
is a net vertical motion of the fractal set of curves of A 
downwards in Fig. 1. We compute the dimension of A to 
be Do = 1.79. 

Arrival times- The first fractal signature can be in- 
ferred from the following argument: Let us fix a point 
P = To in the configuration space. Consider the product 
space of the configuration space and the time axis. In 
this space-time representation, the point P corresponds 
to a 1-dimensional straight line (ro,t), t £ M. The time 
coordinates of the set formed by the intersection of this 
line with the (Dq + l)-dimcnsional time-extension of A 
correspond to the time instants of arrivals of aerosols 
at P. This intersection has a dimension [7] given by 
do ~ 1 + (Do + 1) — {N -1-1), since the time-extended con- 
figuration space is {N + l)-dimensional. Therefore, the 
time series of the instants of arrival of advected aerosols 
at the point P has a fractal dimension do given by (1). 

For practical purposes, we substitute the point P by a 
small codimension 1 object, placed horizontally, which is 
our detector. In the two-dimensional fluid flow that we 
analyze here, such object is a segment of length S ^ L. 
Seen at the scale of the detector, the dynamics of the 



advected aerosols under gravity roughly corresponds to 
a Cantor set of curves moving with some (both time and 
space dependent) horizontal velocity component Vh- The 
fractality of the time series of the instants of arrival of 
aerosols at the detector is measurable at scales r > S/vh- 
To validate numerically this fractal signature, we let 
10^ aerosols, distributed on the attractor within the unit 
square at time zero, evolve until they reach the level 
y = yo = 0. Their instants of arrival as well as their 
X coordinates are then recorded. We then choose a point 
xq and analyze the time series of instants of arrival at 
the segment x G {xo,xo + S), y = 0. The dimension 
dts of this time series should, in face of our theoretical 
findings, be equal to do = 0.79. We find excellent agree- 
ment with this value. For instance, taking 5 = 10~^, for 
xo = 0.2, xo = 0.5 and xq = 0.8, we find, respectively, 
dts = 0.80, dts = 0.77 and dts = 0.80. The numbers 
of points in each of these choices of a;o are, respectively, 
185, 663 and 562. To calculate these dimensions, we have 
used the method described in [12]. This method is much 
more efficient than the usual box counting and is spe- 
cially suited for computing the dimension of subsets of 
the real line. We explain it briefly here: Let us say that 
we want to measure the dimension d of the set M. For 
each point bi € M, let ni(l) be the number of points in 
Af that lie within a distance I of bi. The following scal- 
ing can be shown to hold: (l/n(/)) ~ where d is 
the fractal dimension, and the bracket denotes a uniform 
averaging over all elements of the set. The dimensions 
for the different choices of xo were measured with I = t 
in the range 0.002 < t < 0.05, equally spaced on the 
logarithmic scale. 

Density jumps- The second fractal signature is asso- 
ciated with the spatial distribution P of the advected 
aerosols reaching a line s at a certain level y = yo of the 
configuration space. Let ro £ R^~^ be the spatial coor- 
dinate of a point along a line s at the (A^ — l)-dimensional 
plane y = yo- We define P(ro)dr as the probability that 
an aerosol reaches this plane in a dr-neighborhood of the 
point ro over a finite time interval At. P(ro) is propor- 
tional to the frequency of particles falling near Tq, and 
can in principle be determined experimentally. The idea 
of the second signature is that the downward motion with 
local velocity v caused by gravity is equivalent to a pro- 
jection of the whole fractal pattern shown in Figs. 1, 
3(a), and 4 onto a horizontal line. As illustrated in Fig. 
2, the direction of the projection is tangent to the distri- 
bution at a fractal set of points, which causes discontin- 
uous jumps in the projected measure if At is sufficiently 
large. This means that there is a fractal set of points 
where the density of detected particles has discontinu- 
ities. This is explained more rigorously in what follows, 
and we show that the fractal dimension of A can be ob- 
tained from the dimension of the set of points where P 
changes discontinuously. 

The attractor has an S'PP-like measure [13], which is 
absolutely continuous (discontinuous) along the unstable 
(stable) foliation, and the distribution of points on A in- 
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FIG. 2: (color online) Illustration justifying the fractal sig- 
nature on the density jumps. The black curves represent the 
slice z =const of A at a certain time instant t. The vector field 
V defined on A is represented by the blue arrows. Its tangen- 
cies with A are the extremal points (red dots). The line ^ 
(green) is a I-dimensional line which joins all the extremal 
points of the (x,j/) plane. 
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herits this property. The distribution P is proportional 
to the projected natural measure in the box of size dr 
around rg on the level y = yo integrated over time. The 
measure is not continuous in the configuration space, but 
rather concentrated on the filaments of A. As a result, 
the distribution P can have a discontinuous jump where 
the local velocity v at a point of A happens to be tan- 
gent to the corresponding filament. We call the points 
where this happens extremal points [14]. In order to de- 
termine the dimension of such points, let us first con- 
sider a plane (x, y) of the projected attractor (see Fig. 
2). The extremal points of the main filaments can be 
joined by a smooth line ^ (Fig. 2). The attractor's di- 
mension on this plane is Dq + 2 — N (it is Dq and Dq — 1 
ioi N = 2 and = 3, respectively). Thus, the dimen- 
sion of the extremal points on this plane is that of the 
intersection of a 1 and a Dq + 2 — iV-dimensional object, 
which is 1 (£>o + 2 - iV) - 2 = Z^o + 1 - A^. Next, 
observe that in the product space of the x,y plane and 
the time axis this is a Dq + 2 — iV-dimensional set. The 
points of line s (jj = yo,z =const) on which the distri- 
bution P is defined form a plane in this product space. 
Thus the set of points where the density jumps occur 
in the line s at the fixed level y = yo has a dimension 
2+{Do + 2~N)-3 = Do + l-N'^do. 

The second fractal signature refers to a local property 
in the sense that it is related to discontinuities of a dis- 
tribution. It is, therefore, intrinsically more difficult to 
detect in an experiment, due to the fluctuations induced 
by the finite number of observed aerosols. Such fluctua- 
tions tend to obscure the true discontinuities and lead to 
the consideration of points which do not correspond to 
true discontinuities. In principle, it can, however, be de- 
tected for a sufficiently large number of aerosols. To show 
this, consider the part of the attractor which is shown in 
black in Fig. 3(a). We let 1.7 x 10^ aerosols, forming 
that part of the attractor, evolve until they reach the 
line y = yo = 0. We measure the spatial distribution of 
the aerosols here using a sliding-window method, to be 
described now. We consider the segment (0.88, 0.98) x 0, 
which is divided into 10^ equal boxes s^. We count the 
number of aerosols arriving in each s^. We then define 

Vi ~ X]j=i°° '^j ' * ~ 1; ■ ■ ■ 1 10^ ~ 100, the total number of 
arrived aerosols in 100 neighboring boxes. Next, consider 



FIG. 3: (a) Magnification of the snapshot t = 0.85 of A from 
Fig.l. Both the black and the gray points belong to A. The 
black points correspond to the part of the attractor which is 
used to illustrate the second fractal signature, (b) Density 
jumps Ai as a function of the box index i (in units of 10^) 
(gray) along the line s given by y = j/o = 0. The local maxima 
are also shown (black dots). 
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FIG. 4: Branch of A at f = 0.85 (gray and black points) from 
which the spatial distribution of aerosols is computed in the 
level y = yo — 0. The black subset reaches the level y — ai 
the points where the density jumps (Ai > a) occur. 



Ai = Vi+ioo — Ti, i = 1, • • • , 10^ — 200, the jumps between 
two adjacent intervals, each formed by 100 boxes. Figure 
3(b) shows Ai in the whole range. The local maxima are 
also indicated. We define a threshold a for these jumps 
and consider that the discontinuities in the spatial dis- 
tribution occur at the points corresponding to the local 
maxima of which exceed a. The procedure is robust 
for a in a certain range that, at the same time, allows 
the detection of a reasonably large number of disconti- 
nuities and does not indicate false discontinuities (due to 
the fluctuations). In our numerical experiment, the min- 
imum value of a for which no false discontinuity is de- 
tected is 204. In this case, wc find 51 discontinuities. We 
compute the dimension dsd of this set of discontinuities. 
We obtain d^d = 0.82 over the range 0.0004 < / < 0.01 in 
the interval 0.88 < x < 0.98. The subset of A reaching 
the line y = at the points where these discontinuities oc- 
cur can be seen in Fig. 4. For a = 218, we find 47 discon- 
tinuities, and the corresponding dimension is dsd = 0.81 
over the same range. These dimensions arc in very good 
agreement with the expected value do = 0.79. For these 
computations, we use again the method of Ref. [12]. 
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To summarize, we have shown that important infor- 
mation of the advection dynamics of aerosols in fractal 
sets can be obtained from data measured solely at a given 
level. We have uncovered two independent fractal signa- 
tures of chaotic advection under gravity. We have illus- 
trated our fractal signatures using a flow model with peri- 
odic time-dependence, but our findings arc far more gen- 
eral, since flows with random time-dependence (yet spa- 
tially smooth) also have the properties mentioned here. 
In particular, aerosols advected in such flows often ap- 
proach a random attractor characterized by a well defined 
time-independent fractal dimension [4, 15, 16]. Note that 
the weight of the particles does not play a role in the ar- 
gument. Therefore, the fractal signatures derived would 
also apply for the rising dynamics of finite-size particles 
lighter than the ambient fluid (bubbles). In fact, our frac- 
tal signatures are applicable to fractal chaotic attractors 
in general. 

Our findings might be useful in meteorology, in the case 
where the aerosols are raindrops. Although our chaotic 
advection model neglects possibly important features of 
rain precipitation such as the spatial roughness of tur- 
bulent flows and the distribution of raindrop sizes, the 
fractal signature on the arrival times relics only on the 
time invariance of the dimension of the set where the 
aerosols accumulate. Therefore, if either the aforemen- 



tioned mechanism based on the convergence of aerosols 
in random flows to strange attractors [17] or any other 
mechanism leads to the accumulation of raindrops in frac- 
tal sets with a time-invariant dimension [18], then the 
signature on the arrival times may be used to character- 
ize the fractal clustering of the raindrops. This signature 
could be measurable in precipitation data of rain, for 
instance, with the disdrometer which was used in the ex- 
periment reported in [19]. We note that a disdrometer 
was already used in Ref. [20] and even the correlation 
dimension {D2) of a time series was measured and in- 
terpreted as a sign of irregular distribution of drops in 
space. Finally, we mention the experimental results of 
Ref. [19] showing a fractal dimension for the time series 
of arrival times of raindrops at a disdrometer (Fig. 9(a) 
of the cited reference). Both results [19, 20] arc com- 
patible with the accumulation of raindrops in fractal sets 
with a time-invariant dimension. 
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